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1 ) Problems 


No problems occurred this quarter. 

2) Accomplishments 

a) Geologic anomalies seen on the Death Valley thematic mapper (TM) 
scene were located in the field. Similar anomalies were anticipated 
to be present and subsequently found on the TM scene covering an area 
north of Flagstaff, Arizona. The latter site was also field verified 
and samples from both sites were submitted for analysis to the 
laboratory in Phoenix, Arizona. 

b) A comparison of the information content of a digitally enlarged TM 
image and a digitized black-and-white aerial photograph of the same 
area is being carried out. A northern Virginia TM subscene of the 
Dulles-Reston area was digitally enlarged 7X, smoothed with a 7 by 7 
low pass filter, edge enhanced with a 15 by 15 box filter, and 
registered to the digitized photograph. Approximate pixel resolution 
was four meters. The "blocky" appearance of the TM subscene 
suggested that a smaller (perhaps 5 by 5) box filter for edge 
enhancement would produce a more correct image. 

c) The Dulles-Reston subscene was also compared to a simulated SPOT 
image of the same area for spectral content. The analysis indicated 
that the difference in "color" information present in the three 
equivalent TM bands was not very different from the three SPOT bands. 
High correlation between the equivalent bands of the two sensors and 
little, if any, new information added from band-to-band comparisons 
may have been due to the small size of the test area and a narrow 
range of cover types. 
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d) A paper on TM image processing was presented at the ASP Annual 
Meeting in Washington, D.C. The paper covered some of the techniques 
used to 1) reduce the amount of data that needs to be processed and 
analyzed by either using statistical methods or by combining full 
resolution products with spatially compressed products, 2) digitally 
process small sub-scenes to both improve the visual appearance of 
large scale products and to merge different resolution image data, 
and 3) evaluate and compare the information contents of the different 
three band combinations that can be made using TM data. Results from 
the use of techniques for assessing information content indicate that 
for some applications, the added spectral information over MSS is 
more important than the Increased spatial resolution. This is 
because there is "new" information in the TM bands that were not part 
of the MSS system. For example, the potential use of two three-band 
TM combinations for extraction of soil moisture and ground water 
depth Information in arid and semi-arid areas shows promise and 
warrants further research. 

e) The statistical information generated earlier from TM data for the 
Dyersburg and San Francisco areas was checked. The earlier results 
were found to be correct and have been added to the statistical/ 
Optimum Index Factor (OIF) package with other TM image results. The 
San Francisco TM image has also been used to evaluate water 
enhancement methods developed for MSS data. The results were not as 
good as was expected because of the amount of noise (striping plus 
other high frequency noise patterns) that was present in TM band 1. 

In order to improve the results, it is necessary to go back to the 
step before stenciling (removal of non-water pixels) and use a second 
order convolution filtering technique to remove the noise before 
applying the stretch, smoothing, and color coding algorithms. 

f) Analysis of the Washington, D.C. TM image was requested by NMD in 
Reston. The statistical information was needed to generate OIF 
rankings and to perform selective principal component analysis on the 
data. The OIF technique ranked combination TM 1 , 4, and 5 as 
number 1 and TM 1 , 2, and 3 as number 20 (lowest). 

g) Black-and-white and color prints of processed TM data of the Silver 
Bell copper mine area were sent to General Electric (GE). GE 
requested, and NASA approved, use of some of the prints in their 
Landsat-TM brochure prepared for distribution during the March 1 
launch of Landsat D. 

3) Significant Results 

Extensive work has been done using the Optimum Index Factor (OIF) for 
choosing TM band for color compositing or reducing the amount of digital 
data to be processed. See publication attached. 
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4) Publications 

Chavez, P. Jr., Guptill S. C., and Bowell, J. A., Image Processing 
Techniques for Thematic Mapper Data: Proceeding of the ASP Spring 

Conference, Washington, D.C., March 1984. 

5) Recommendations 
None 

6) Data Utility 
None 
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ABSTRACT 

The Thematic Mapper (TM) Imaging system onboard Landsat 4 
provides new and much more comprehensive data than the data 
collected during the past 10 years by the Landsat multl- 
spectral scanner (MSS) system. A major challenge In working 
with the TM data Is to process and analyze In an efficient 
manner the much larger volume of data resulting from the 
higher spatial resolution and the Increased number of spec- 
tral bands. Processing techniques have been used on TM data 
over several sites to (1) reduce the amount of data that 
needs to be processed and analyzed by using statistical 
methods or by combining full-resolution products with spa- 
tially compressed products, (2) digitally process small sub- 
areas to l^rove the visual appearance of large-scale pro- 
ducts or to merge different-resolution Image data, and 
(3) evaluate and compare the Information content of the 
different three-band combinations that can be made using the 
TM data. Results Indicate that for some applications the 
added spectral Information over MSS is even more Important 
than the TM*s Increased spatial resolution. 


INTRODUCTION 

The Thematic Mapper (TM) Imaging system onboard Landsat 
which was launched on July 16, 1982. provided new and much 
more comprehensive data than Che data collected during the 
past 10 years by the Landsat multlspectral scanner (MSS) 
system. The MSS system collects data In four spectral bands 
at a pixel size of approximately 80 meters, and the TM sys- 
tem collects data In six spectral bands having a pixel size 
of approximately 30 meters and one band (thermal) with a 
pixel size of about 120 meters (see Table >1 for a compar- 
ison between MSS and TM characteristics). A major chal- 
lenge In working with the TM data over the MSS Is to pro- 
cess and analyze In an efficient manner Che much larger 
volume of data resulting from the Increased spatial and 
spectral resolution. 


HEASURg 

TABLE 1 

LAHOSAT-4 INSTRUMENT COMPARISON 
MULTISPECTRAL SCANNER fMSSl 

THEMATIC mapper fTH) 

KOMIMAL GROUNO 

60 s 80 HETgR 

10 X 30 METER 

RESOLUTION 

SUAIH WIDTH 

186 KM 

18S KM 

39ECTXU. BANOS (W1) 

— 

0.4S TO 0.52 


0.5 TO 0.6 

0.52 TO 0.60 


0.6 TO 0.; 

0.63 TO 0.69 


0.7 TO 0.8 

0.76 TO 0.90 


0.8 TO 1.1 

1.55 TO 1.75 


— 

2.08 TO 2.35 



10. « TO 12.6 

A/0 CONVERSION (BITS 

6 

8 

PER PIXEL) 

APPROXIMATE NUMBER OP 

S2 HILUON 

300 MILLION 

PIXELS PER SCENE 

(ALL BANOS) 



DATA RATE 

IS MBPS 

85 MBPS 


(TABLE FBOH 0. H. SMITH. 1982) 


A variety of processing techniques have been used on TH 
data over several different sites to address three topics. 
The first topic is to reduce the amount of data that needs 
to be processed and analyzed by either using statistical 
nethods or by combining full-resolution products with spa- 
tially compressed products. A second topic is to use digi- 
tal processing on small sub-areas to improve the visual 
appearance of large scale products or to merge different- 
resolution image data. The final topic is to evaluate and 
compare the information content of the different three-band 
combinations that can be made using the TM data (excluding 
the thermal band) . 

The digital TM data used in this study was obtained from 
NASA's Goddard Space Flight Center through the EROS Data 
Center under two separate NASA TM investigations. The data 
was in P-tape format, with some radiometric and geometric 
corrections performed at Goddard (Salomonson and others. 
1980) . This paper does not deal with the geometric accur- 
acy of the TM images and the interested reader is referred 
to papers by Blanchard and Weinstein (1980) and Wrigley and 
others (1983) for discussions on this matter. 

Data Processing and Analysis 

One of the main products generated from the TM data is a 
color composite using three of the six bands (excluding the 
thermal band) . Because a total of 20 combinations can be 
made from the six TM bands taken three at a time (ignoring 
the permutations of bands and display colors), deciding 
which combination contains the most information on the basis 
of visual analysis can be difficult and time-consuming. The 
problem is more difficult with the 15 different ratios that 
can be generated from the six TM bands because they can be 
combined in 455 different ways (15 ratios taken three at a 
time). A technique called the Optimum Index Factor (OIF) 
developed earlier by one of the authors was used to rank the 
20 possible combinations in the test sites studied (Chavez 
and others. 1982). The OIF technique ranks all the combin- 
ations possible based on the amount of correlation and the 
total variance present between the various data sets being 
used (e.g.. six TM bands). The algorithm used to compute 
the OIF value for any subset of three bands is as follows: 


3 3 I 

OIF = E SDl / £ CCJ 
l»l Jail 

where. SOL = scandard devLacion for band 1 

I CCJ I » absolute value of Che correlaclon coefflclenc 
between any two of the three bands being 
used (Chavez and other, 1982). 

The combination having the largest OIF value should have Che 
moat information (as measured by variance) with the least 
amount of duplication. Often, the combinations chat are 
within 2 to 3 ranking positions of each ocher appear similar 
and distinct feature differences cannot be detected. This 
is because there is little difference between closely ranked 
combinations as far as total information content is con- 
cerned. Table 2 shows the correlation matrix and the OIF 
rankings for several sites using the TM data. 

The information shown in Table 2 was derived from TM data 
analyzed from a wide range of different test sites. The 
Death Valley data represents an arid environment; the 
Washington, O.C. area is an agricultural and urban setting; 
Che San Francisco site includes agricultural, highly urban, 
and water settings; and finally, the Dyersburg, Tennessee 
site is mostly an agricultural setting. From the analysis 
of Che data and Che OIF rankings it was found that, in gen- 
eral, Che three-band combinations chat included one of the 
visible bands (TM 1, 2, or 3) and one of Che longer wave- 
length infared bands (TM S or 7) with TM 4 were usually 
ranked high by the OIF technique (i.e., having the most in- 
formation with the least amount of duplication). This hap- 
pens because of the high correlation chat exists between TM 
bands 1, 2, and 3 and between TM bands 5 and 7 (see 
Table 2). 

Additionally, in an attempt to maximize the amount of in- 
formation contained in a final chree-band/component com- 
bination and to minimize the information lost due to not 
using Che ocher three bands, selective principal component 
analysis can be used (Chavez and others, 1982). Selective 
principal-component analysis involves using only highly 
correlated subsets or pairs of bands as input to principal- 
component analysis and not using all the bands simultane- 
ously. Because the bands selected for principal-component 
analysis are so highly correlated most of the information 
is mapped into the first component. With the TM data pro- 
cessed so far, Che subsets or pairs of bands used as input 
to selective principal-component analysis have been TM bands 
1, 2. and 3 as one group and TM 3 and 7 as the second group. 
The first component of each of the two results is then com- 
bined with TM q to make a new data sec that can be used as 
input for further processing (for example, contrast and edge 
enhancement for color compositing or input to digital 
classification) . 

In areas where this has been done between 96 to 98 percent 
of Che information (variance) contained in the six bands 
was mapped to the new three-band/component combination. 
Silver Bell, Arizona was one of the sices processed by 
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selective pclnclpal-componenc analysis and the results are 
shown In Table 3. The user does have to be careful to checlc 


TABt^ : 


TH SIX-BAMO CQR«EL*nON KATAIX FOR THE SILVER BELL AREA 
(IHt • THERMAL BAND NOT USES) 
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SEXJECTEO SUBSET OE BANDS FDR 'SELECTIVE* ?R1NCI?AL 
COKPOHENT ANALTSXS BASED ON HIGH CORRELATION 

I. Tin PCI3 96.81 Of VARIANCE 

TM PCM 1. 61 OP VARIANCE 

TM PC33 — ^ 0.61 OP VARIANCE 

II. TM5 PC12 " ■ ' 98. *1 OP VARIANCE 

TH7 PCM ■ ' ' 1.61 OP VARIANCE 

III. USE PCI3 . PC13 . and as THREE-COMPONENT SATA BASS 
FOR ANALXSES. CLASSIFICATION. OR COLOR COMPOSITING. 
CONTAINS A HINIMUM OP 98.31 OP SIX BAND VARIANCE. 



the components/data not used to oalce sure that Information 
of Interest has not been Isolated or mapped to one of the 
lover components. This did occur In the Silver Bell data 
where most of the Information about an alteration sone of 


geologic Interest was mapped to 
S and 7 ( see f Igures 1 and 2 ) . 



component number two of TM 
(This component also 


Figure la . Landsat <* TM 
Image of band 5 (1.55 to 
1.75 m) with a linear con- 
trast stretch of 9-108 (99 
ON range). The Silver Bell 
copper mine Is located In 
the center of the Image. 


happened to be highly correlated with the ratio of TM 5 to 
TM 7 (correlation coefficient of -.93)). Therefore, If the 
user Is Interested In one specific cover type Instead of 
maximizing the information content over the entire area, 
selective principal-component analysis must be used with 
care, as should regular principal-component analysis. For 
information on principal-component analysis and other exam- 
ples of Its use, Che reader is referred to Blodgec and 
others (1978) and Pirkle and others (1980). 



Figure lb . Landsat 4 TM 
Image o£ band 7 (2.08 Co 
2.35 m) with a linear 
contrast stretch of 4-57 
(53 ON range). 


Another method that can be used to reduce Che amount of data 
Chat needs to be processed and analyzed is to use a combin- 
ation of full- and compressed-resolution data sets. For 
many applications that require high-frequency or fine- 
detail Information only one TM band In black and white with 
an edge enhancement Is needed (for example, mapping of lin- 
ear features). There are some applications where the added 
spectral Information of TM over Che MSS (7 bands vs 4 bands) 
Is much more Important chan the Increased spatial resolution 
(30 metrers vs 30 meters). In a test conducted using Che 
Death Valley TM data where Che pixel size was compressed to 
60 and 120 meters, there are color anomalies (discussed 
later In this paper) that were still visible on Che 120- 
meter pixel size Image made from TM Bands 1. 4. 5 but were 
not visible on Che composite made from TM bands 1, 2. and 3 
(natural color — NC) or TM bands 2, 3. and 4 (similar to 
MSS) at full resolution. In this case, the compressed data 
were complemented by an edge-enhanced black-and-white Image 
of TM 4 at full resolution for fine-detail mapping. 

Research Is continuing Into various ways to compress the 
digital data. Methods being compared Include dropping 
pixels, averaging pixels, and averaging only where the first 
difference within the area to be compressed Is small or 
using Che pixel value within Che area where Che largest 
first difference occurs If It exceeds a sec threshold. By 
doing this we hope to be able to reduce Che Image size but 
still retain full resolution In areas where It Is needed. 

AC times, rather chan being concerned about Che large volume 
of data chat exists for processing, a user wants Co take 
data for a small area and enlarge It or Integrate It with 
higher resolution data. Digital enlargement combined with 
spatial filtering can be used to accomplish this. A digital 
Image (MSS, TM. or digitized photo) that Is optically en- 
larged close Co or past Che limits supported by Its pixel 
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Figure 2a . Principal com- 
ponenc number 1 of TM 5 and 
7 with a linear contrast 
stretch of 9-120 (111 ON 
range) . The amount of the 
total variance of TM 5 and 
7 mapped to principal com- 
ponent number 1 was 98.36 
percent. 




Figure 2b . Principal com- 
ponent number 2 of TM S and 
7 with a linear contrast 
stretch of 69-89 (20 ON 
range). The amount of the 
total variance of TM S and 
7 mapped to principal com- 
ponent number 2 was 1 . 64 
percent. This image also 
shows a noise pattern con- 
tained In the TM data. The 
areas showing up In the 
darlcer tones correspond 
closely to an alteration 
zone previously mapped by 
the Geosat committee and JPU 
using TM simulator data. 



Figure 2c . This shows the 
results of the TM 5 to TM 7 
ratio with a linear contrast 
stretch. Notice the corre- 
lation between It and prin- 
cipal component number 2 
(correlation coefficient of 
-.89) . 







Figure 2d . Sfune Image shown 
In figure 2b (l.e., princi- 
pal componenc number 2 of 
TM 5 and 7) but with a small 
smooching filter to suppress 
the noise and a harder con- 
trast stretch to show the 
alteration zone area better. 


Figure 2e . Same Image shown 
in figure 2c (l.e., TM 5 to 
TM 7 ratio) but with a small 
smooching filter to suppress 
Che noise and to compare 
with principal componenc 
number 2 shown above. This 
ratio and principal compo- 
nenc number 2 above have a 
correlation coefficient of 
-.93. 




cesoluClon will have a "bloclcy” appearance and often a fine 
screen pattern (caused by an anomaly In the hard-copy film 
writer) can be seen (see figure 3a). To solve this problem 



Figure 3a . Portion of TM band 5 Image showing the 
Washington. D.C. Dulles airport area. This Image represents 
the original data without any digital enlargements. The 
bloclcy appearance and screen pattern discussed In the text 
can be seen In this product, especially with a 8X lens. 


we have used a method developed by one of the authors In his 
Investigation of geometric resampling. The process Involves 
the combination of digital enlargement by pixel duplication 
followed by spatial filtering for smoothing and. If neces- 
sary, edge enhancement. The main reason that an Image ap- 
pears bloclcy and the screen pattern becomes visible Is that 
when an Image Is optically enlarged to a scale that Is close 
to or past the limits supported by Its pixel resolution, the 
Individual pixels become so large that they no longer In- 
teract visually to create patterns and features. Each pixel 
starts to be viewed separately and Is seen as an Individual 
feature or pattern. To minimize this problem the number of 
pixels In the Image Is Increased so that Che- less optical 
enlargement Is needed to get the Image to the desired scale. 


resulclng Ln a smaller final pixel size. However, If only 
pixel duplication Is used the Image will still look block/ 
at Che desired scale (see figure 3b). The block/ appearance 



Figure 3b . Same Image shown In figure 3a but with a 3X 
digital enlargement. The enlargement was done using simple 
pixel duplication so that one pixel became a 3 by 3 array 
where all nine pixels were assigned Che same digital value. 


may be more obvious Chan In the original Image because one 
pixel becomes an NxN pixel array with the same brightness 
(for example. 3x5 In Che case of a 5X digital enlargement). 

To eliminate this effect a smoothing filter needs to be 
applied to Che data. The size of Che filter Is critical 
(see figure 3c). and must be equal to the size of the 
digital enlargement In order to not over- or under-smooth 
the data. Also, by using a smoothing filter that Is Che 
same size as the digital enlargement, the values of Che 
pixels at Che center of the enlarged arrays will not change 
(for example, pixel at middle of 3x5 array using a 3X 
enlargement). This has the same effect as separating Che 
original pixels by Che digital enlargement factor and then 
using a two-dimensional smoothlng/lnterpolaclon technique 
to fill In Che missing holes (l.e., resampling). 

Digital enlargement by pixel duplication can be considered 
a form of nearest-neighbor resampling when applied to an 
Image chat already has Che desired geometric character- 
istics. "Die two-dimensional smooching filter Is a form of 
Interpolation similar to bilinear or cubic convolution. 




Figure 3c . Sane Image shown In figure 3b but with a 3 by 3 
smoothing filter applied to the data. It Is Important that 
the size of the smoothing filter be equal to the size of the 
digital enlargement In order to not over or under smooth the 
data. 

This method Is more straightforward and computationally 
efficient than resampling the original Image using cubic 
convolution to a higher resolution. Also, the procedure 
can be repeated with different enlargement factors If the 
user desires several different scales or sub-areas from the 
same Image. 

This enlargement technique has been used on both Landsat 
MSS (5X) and TM (3X and 7X) data with very encouraging re- 
sults. The MSS data held up well at a scale of 1:100,000 
and the TM at a scale of 1:50,000. The method also helped 
suppress some of the random noise effects often seen In en- 
larged digital Images and produced more homogeneous patterns 
(this can be seen by color coding the before and after 
Images using the same color table). The method expands the 
pixel spacing and Interpolates a smooth transition between 
the brightness values of two pixels Instead of having an 
abrupt change. This Is similar to what bilinear, cubic 
convolution, or restoration does when resampling MSS data 
to 50-meter pixels. In the 3X enlargement of TM data the 
resultant Image looked slightly defocused. To reduce this 
apparent defocuslng algorithm an edge enhancement was ap- 
plied to the data using a kernel size that was approxi- 
mately twice the digital enlargement (that Is, digital 


enlatgemenc of 3x3 with an edge enhancemenc Vernel of 7X7, 
see figure 3d). In Che MSS and TM examples which used a 5X 
and 7X enlargement followed by 5x5 and 7x7 smooching 



Figure 3d . Same Image shown In figure 3C but with a 7 by 7 
edge enhancemenc. When the digital enlargement Is small 
(e.g., 3X) Che resultant Image can appear slightly out of 
focus so that an edge enhancemenc will be needed to 
sharpen-up the Image. 


filters. Che resultant Images did not look as defocused. 

This may be due to the Increased number of new pixels 
placed among the original pixels thus reducing the 
Incremental brightness transition between pixels. The 
products generated by this method will be compared with 
some made by directly digitally enlarging an Image using 
cubic convolution resampling. 

The Image analyses comparing the TM bands which are similar 
to Che MSS bands (TM 2, 3, and with other TM band com- 
binations. which Include some of the new additional spec- 
tral bands. Indicate that there Is new or more information 
In the TM data set. To evaluate only the spectral differ- 
ences and not be Influenced by the difference In spatial 
resolution. Che comparisons have been made using only the 
TM and not the MSS data (this will be done In future work). 
TM bands 1. 2. and 3. which approximate the visible part of 
Che spectrum and generates a composite that has colors close 



CO nacural color, were also compared with the TM 2, 3. and 4 
combination. Table 2 shows that both the natural-color and 
MSS-equlvalent band combinations are usually ranked quite 
low by the Olf method. 

Visual analyses of several TM band combinations In the 
Death Valley Image have produced some Interesting results. 
There are some surflclal anomalies that are readily apparent 
In this area on TM band combinations 1. 5. 7 and 1, 5. 

combinations ranked 1 and 3 by the OIF (see Table 2). On TM 
band combinations 1. 2. 3 and 2. 3, 4 (natural-color and 
NSS-equlvalenc) . there are no signs of these same anomalies. 
Several of these anomalies were field-checked and once 
Identified. It was correctly predicted that similar ano- 
malies would appear In the same TM combinations In an Image 
covering an area Just north of Flagstaff, Arlz. These areas 
have also been field-checked, and the results of that In- 
vestigation are currently being prepared for publication. 

Further analyses of some Death Valley Images have Identified 
another potential use of the TM band combinations. On some 
TM band combinations It appears that there may be a rela- 
tionship between soil moisture/ground water depth and the 
colors seen on the resultant Image product. In the desert 
Death Valley Image TM band combinations that Include either 
TM bands S or 7 or both portrayed as green and red. respec- 
tively, In^color composite, seem to map areas Influenced by 
high soil moisture or shallow ground water In various shades 
of blue. Several playas In the Image appeared as light- to 
medium-blue while others remained white. Also, drainage 
areas along the playas were usually darker than their sur- 
rounding light-blue areas. Indicating lower brightness 
levels In TM bands 5 and 7. Snow on the mountains sur- 
rounding the Death Valley area appeared In the same Image 
products as very dark blue, indicating again low brightness 
values In TM bands 5 and 7 . 

One possible explanation Is that TM 1, a visible band not 
affected by water absorption, was exposed through the blue 
filter and TM i and TM 7 were exposed through the green and 
red filters, respectively. Because TM 5 and TM 7 are af- 
fected by water absorption, areas of higher soil moisture 
will have lower digital values than their surrounding areas 
with lower soil moisture. Areas appear In shades of blue 
depending on the amount of soil moisture because the 
brightness values are less In TM 5 and 7 (used as green and 
red In the color composite) while they stay relatively con- 
stant with the surrounding areas In TM 1 (used as blue In 
the composite and not affected by soil moisture). This 
color difference can be seen less clearly In composites made 
with TM 1. 2. 5 and TM 1. 2. 7. 

It must be pointed out that these soil moisture correlations 
are only preliminary, and no field calibrations have been 
possible due to data delivery problems encountered with the 
TM system. Most Investigators are working with data col- 
lected during time frames that could not be predicted. The 
current observations do warrant further Investigation when 
the new TM system Is launched, offering opportunities for 
appropriate field calibration. 


SUMMARY 


The techniques of OIF ranking of band combinations. Image 
enhancement, and spatial compression have been used to pro- 
cess and analyze efficiently large volumes of TM data. From 
the digitally compressed data and other analyses. It seems 
that there ace some applications vhere the added spectral 
Information In Che TM data Is more Important than the In- 
creased spatial resolution, as compared to Che MSS. Also, 
digital enlargement and smooching algorithms have been used 
successfully to generate large-scale TM and MSS products and 
to merge data sets with different resolutions. 

The Image analyses done under this project Indicate that 
there Is new Information In Che added TM bands Chat was not 
provided by the MSS system. The potential use of TM data 
to extract Information about soil moisture and ground-water 
depth warrants further research and future data collection 
where field calibration Is possible. 
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